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Despite many years of study of glass-forming liquids, 
the most suitable paradigm in which to discuss the 'glass 
transition' and its associated dynamical phenomena re- 
mains controversial. In recent years, there has been 
progress 0, 0, 0J, driven by the idea that the dynami- 
cal properties of glass-formers may be characterised by 
a zero temperature dynamical fixed point 0]. This is in 
contrast to the predictions of other theories which involve 
a finite temperature singularity in the dynamics 0, Q or 
the thermodynamics of the relevant system 0, Glial- 

However, there remains an important qualitative dif- 
ference between physical glass-formers and the mo dels, 
such as the Fredrickson- Andersen (FA) [91 or Eastjlfjj 
models (see 0| f° r a review) studied in [l| 0, 0: the 
crystal phase is completely absent from these models. As 
a result, they are relevant if the glass- forming liquid is in 
a long-lived (but metastable) supercooled phase. Fur- 
ther, interplay between the thermodynamic singularity 
associated with the transition to the crystal and the dy- 
namical fixed point associated with the glass is certainly 
possible, but the FA and East models cannot capture 
these phenomena since their thermodynamics are triv- 
ial. There has been recent work[l2l fl3| (see also Q) 
investigating these issues with regard to glassy systems, 
although without reference to a glassy fixed point at zero 
temperature. 

The extent to which the FA and East models can be 
viewed as pictures of real glasses is therefore contingent 
on two main assumptions. Firstly, the behaviour of the 
supercooled state should not be affected by the proximity 
of the freezing transition, since the FA and East models 
regard glassy slowing down as a purely dynamical effect, 
not requiring a thermodynamic transition. Secondly, the 
'mobility field' represented by the spins in these models 
should emerge naturally from atomistic degrees of free- 
dom of the glass-former. 

One class of models in which this latter effect is 
demonstrated are the two-dimensional plaquette models 
flH ITfjl IItI , in which an effective dynamical constraint 
emerges naturally from a simple spin model. There are 
free excitations in the spin field which are naturally inter- 
preted as a mobility field. At temperatures lower than 



the glassy onset temperature, T a 0, 0, the dynamics 
are strongly heterogeneous, and slow down rapidly with 
decreasing temperature. 

In this paper, we address the other issue mentioned 
above: how are dynamics of metastable supercooled 
states affected by the presence of the freezing transition? 
We study dynamics in the eight-vertex model, whose 
thermodynamics were solved by Baxter [l^. The model 
may be treated as a generalisation of the square plaque- 
tte model 0, 0| ; the effect of this generalisation is to 
to introduce a (second order) phase transition to an or- 
dered state at a finite temperature T c . We identify this 
with a freezing transition, and investigate the dynam- 
ics around this transition. The transition temperature 
T c may be varied with respect to the glassy onset tem- 
perature T Q . This fact, together with the exactly solved 
thermodynamics, gives an extra degree of control to our 
simulations. 

We will show that we may prepare long-lived 'super- 
cooled' states below T c , whose dynamics are controlled by 
the effective dynamical constraint of the plaquette model, 
and are not affected by the freezing transition for times 
shorter than the lifetime of the supercooled state. 

To be more precise, the dynamics of the system within 
the supercooled state resemble those of strong glasses, 
and arise from diffusing point-like excitations in a mobil- 
ity field 0,0, El- Considering supercooled states at differ- 
ent temperatures, we find that they obey dynamical scal- 
ing consistent with a zero temperature fixed point. The 
presence of the ferromagnetic transition means that these 
states have a finite lifetime, but it does not affect the dy- 
namics on timescales shorter than this lifetime. This is 
consistent with the assumptions made when modelling 
glass-formers using models without an ordering transi- 
tionfjiai. 

The form of the paper is as follows: in section^ we de- 
scribe the model, identify relevant energy scales and their 
hierarchy, and discuss the relation between the spins of 
our model and the atomistic degrees of freedom of a phys- 
ical glass former. We discuss the nature of the ordered 
and disordered phases of the model system in section ITT1 
we then use this information to interpret simulations of 
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the dynamics of the model in section IIIII Finally, we 
summarise our results in section II VI and discuss their 
significance for models of the glass transition. 

I. THE MODEL 

The zero-field eight vertex model, solved by Baxter[l^ 
may be expressed either in terms of its original vertices, 
or as an Ising model with Hamiltonian: 

H = ^2ij [ — D<Tij<Ji+l.j(TiJ + l<Ti+lj + l 

— J (oyOi+ij+i + <7i+l,j<7i,j+l)] (1) 

where the {(7^} are Ising spins on a square lattice. 
We note that the Ising coupling is between next near- 
est neighbours on the square lattice: at D = there 
are two independent sublattices, with Ising coupling, J 
within each sublattice. There is a transition at J = 
(T/2) sinh _1 (l) to a fourfold degenerate ordered state 
(there are two ferromagnetic and two antiferromagnetic 
ground states, related by flipping all the spins on either 
sublattice). As D is increased from zero, the lattices 
become coupled, and the transition moves to a higher 
temperature: the critical temperature T c satisfies: 

sinh(2J/T c ) = exp(-2D/T c ) (2) 

The transition to an ordered state occurs for all finite J: 
we also note that if D > J then the transition tempera- 
ture will be much larger than J. 

Thus far we have considered only static (thermody- 
namic) properties of the eight-vertex model. In order to 
study the time evolution of the this model, we must spec- 
ify dynamical rules. We use simple spin flips with rates 
given by Glauber dynamics. We refer to the combina- 
tion of the Hamiltonian and the dynamical rules as the 
spin-flip eight vertex (SEV) model. 

If we set J = in the SEV model, we arrive at the 
(two dimensional) plaquette model0, |2(], Hl|. In this 
limit there is no ordering at any finite temperature: all 
two point static correlations vanish. This may be most 
easily demonstrated by noting that if J — then the 
Hamiltonian is invariant under flipping all of the spins 
in any row or column of the square lattice. The dynam- 
ics of this model are dominated by a zero temperature 
dynamical fixed point for temperatures T that are small 
compared to D. Since we are studying slow dynamics we 
work throughout at T < D. 

We have now identified two temperature scales in 
the problem: the onset of glassy dynamics occurs at 
T = T Q D and the critical point in the system is at 
T = T c . If T c > T a then we expect the slow dynam- 
ics to be observable only in the ordered phase: the more 
interesting case is T Q > T c , in which case the dynamics 
are slow near the transition, and we may investigate the 
effect of the effective kinetic constraint as we cool the 
system through T c . We therefore work at T c < D: from 
(J2J), this means that J < T c . As a result we have the 



hierarchy J < (T,T C ) < D which is obeyed throughout 
this paper. 



A. Relation of this work to physical glass formers 

Before investigating the SEV model more closely, we 
establish the relationship between this model and phys- 
ical glass-formers: it is not obvious at first sight pre- 
cisely how a model of Ising spins should be related to a 
atomistic system. The spin variables represent the mi- 
croscopic degrees of freedom of the glass-former. This 
is distinct from the more heuristic approach taken in 
the FA and East models in which the spins represent 
a coarse-grained 'mobility field'. The SEV model is more 
similar to the plaquette model0] in that the effective ki- 
netic constraint responsible for the critical slowing down 
at zero temperature is not inserted explicitly, but arises 
from the combination of the Hamiltonian and simple 
spin- flip dynamics. In section lTTTl we will comment briefly 
on how the dynamics may be interpreted within a 'mo- 
bility field' picture like that of the FA model. 

In the previous section we identified the two tempera- 
ture scales in the model as the glassy onset temperature 
T D ~ D and the transition temperature T c < T Q . These 
separate the behaviour of the system into three regimes. 
We argue that the high temperature phase of the SEV 
model with T > (T ,T C ) resembles a liquid-like state of 
the atomistic system, since it lacks any two point corre- 
lations between the spins. 

The second regime is T c < T < T Q : there are still no 
static correlations between the spins but there are strong 
dynamical correlations. This state resembles a viscous 
liquid whose relaxation time is large compared to micro- 
scopic timescales. We emphasise that the crossover be- 
tween this regime and the high temperature behaviour is 
smooth: there is no sharp transition at T a . In the viscous 
liquid, the atomistic degrees of freedom are 'jammed' over 
large regions of the system: relaxation in these regions 
is frustrated by large energy barriers. Dynamical het- 
erogeneity then arises naturally, due to the presence of 
mobile regions where the energy barriers are smaller than 
average. There are very many paramagnetic states in the 
spin system, even at low temperatures (compared to D): 
these resemble the many possible jammed states of the 
glass former. 

Having argued that the paramagnetic phase of our 
model is liquid-like, and shows glassy behaviour at low 
temperatures, it is natural to interpret the transition in 
the model as a a freezing transition. We identify the 
ferromagnetic phase with the crystalline states of the 
glass-former. As the temperature is lowered through T c 
the entropy falls rapidly as the very many paramagnetic 
(jammed) states are now thermodynamically unstable 
with respect to the ferromagnet. The effect of the dy- 
namical fixed point on the transition between paramag- 
net and ferromagnet is the main subject of the following 
sections. In particular we show that 'supercooling' of the 
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FIG. 1: Sketch showing the relation between spin states, ex- 
cited plaquettes and domain walls. Up (down) spins are rep- 
resented by filled (empty) circles. Domain walls form closed 
loops with excited plaquettes (energy cost 2D) at each iso- 
lated corner (there is no energy cost associated with wall 
crossings). A domain wall costs energy 2.1 per unit length 
when running though a ferromagnetically ordered region. In- 
teractions between domain walls reduce this tension in disor- 
dered states. 



paramagnetic state is possible as long as T c <g; D. 



II. STATIC PROPERTIES OF THE SEV MODEL 

We now discuss the microscopic structure of the ferro- 
magnetic and paramagnetic states in the model of equa- 
tion Q). This will allow us to identify the relevant cor- 
relation functions for our study of the dynamics of the 
system. We describe the paramagnetic state in terms of 
small deviations from the behaviour of the model with 
J — 0, and the ferromagnetic state in terms of excita- 
tions in an ordered background. This will lead us to 
interpret the transition in terms of free defects above T c 
that become confined at the transition, forming compos- 
ite excitations. We will also show that these descriptions 
are valid even rather close to the transition, despite be- 
ing based on expansions around the fully ordered or fully 
disordered states. In other words, the critical region is 
very narrow. 

We begin by recalling some results for the J = 
limit of the model [20L Elj . At J = 0, we write pi — 
<Jij<Ji.j+i<Ji+i,jtJi+i,j+i, and the Hamiltonian reduces to 



H 



j=o 



Pi 



(3) 



where the Pi are Ising-like degrees of freedom, defined 
on the plaquettes of the square lattice. In the thermo- 
dynamic limit, these plaquettes are independent degrees 
of freedom, that define the state of the spin system, up 
to transformations that flip all the spins in any row or 
column (leaving Hj = q invariant). We observe that this 
results in a ground state entropy proportional to the lin- 
ear size of the system, L (there are N = L 2 spins in the 
system). 

In finite systems, the presence of boundary conditions 
may impose constraints on the plaquettes. For example, 




FIG. 2: Typical states of the spin system above T c (left) and 
below T c (right). Some of the excited plaquettes are identified 
by arrows. They are 'free' above T c , but 'confined' on corners 
of rectangular excitations in the ferromagnetic phase. 



imposing periodic boundaries on the spins means that 
the number of excited plaquettes in all rows and columns 
must be even (we believe that this fact led to the strong 
finite size effects seen in pH). 

As discussed in |2l| , the low temperature states of the 
model with J = are best interpreted in terms of closed 
loops with excited plaquettes at each vertex (see figure^). 
If we move across the lattice, any spin flip is accompanied 
by our crossing the perimeter of a loop. If J = then 
each plaquette is independent: the vertices of the loops 



are a free lattice gas with density (1 
free energy per site is simply 



J1D 



/Ty 



/j = o = -Tlog[2cosh(D/T)] 



The 



(4) 



At finite J, we make use of Baxter's solution of the 
eight vertex model^^- ^ n appendixEI we show that the 
free energy per site for T 3> T c is given approximately by 



m = -Tlog \e D ' T cosh(2 JIT) + e- D ' T 



(5) 



Since we work exclusively at J < T, equation JSJ) is rather 
close to the J = expression 

We focus on two correlation functions, the density of 
excited plaquettes n p , and the density of broken Ising 
bonds, rib- Their definitions are: 



-(1 - (o'ij'0'ij+lO'i+ljCri+lj + l)) 



(6) 



n b = j(2 - ((7yCr. i+ lj + l + Vi+lj&i.j + l)) (7) 

In the representation of figure ^ the concentration of 
vertices is given by n p . The parameter nb is related to 
the total perimeter of the closed loops of that figure, and 
measures the spatial ordering of the vertices. The free 
plaquettes observed at J — have Ub = (1/2). As nb 
is reduced, the reduction in the loop perimeter starts to 
constrain the positions of the vertices, and spatial corre- 
lations appear. 

The internal energy per site is given by 



(H/N) = 2Dn p + AJn b - (D + 2 J) 



(8) 
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so we have ^ = {X + df/dD)/2 and n b = (2 + df/dJ)/4. 
Using J5J), the paramagnetic state has 



rib 



[l + e 2D / T cosh(2J/T)}- 1 
(1/2) - Smh(2J/T) 



cosh(2 J/T) + e - 2D / T 



(9) 
(10) 



We see that introducing J leads to small negative cor- 
rections to the J = values of n p and n&. However, 
for T < D, the concentration of excited plaquettes is 
still approximately c = e ~ 2D / T and these plaquettes are 
only weakly interacting since n& ~ (1/2). A typical con- 
figuration is shown in figure [2 there are no two point 
correlations between the spins, but the loop vertices are 
dilute since T < D. 

We now consider the system at temperatures lower 
than T c . As shown in appendix E] a good approximation 
to the free energy per site in the ferromagnetic phase is 



AO) 

Jfm 



-D-2J 



Te -8D/T 



16 sinh 2 (2 J/T) 
which is valid for J < T < T c < D. In this regime 

e -8D/T 



4 sinh 2 (2 J )T) 



,-SD/T 



rib 



16 sinh 2 (2 J )T) tanh(2 J/T) 



(11) 

(12) 
(13) 



For excitations in a ferromagnetic background, nb is di- 
rectly related to the perimeter of the closed loops shown 
in figure ^ Equations 112|) and (|13|l are therefore con- 
sistent with rectangular excitation loops with an excited 
plaquette at each corner (total energy cost 81?) . The ex- 
pectation of the loop perimeter is approximately (T/2J). 
This is much smaller than the typical spacing between 
loops, given by (8J/T)e 4D / T . This situation is sketched 
in figure 

The density of loops is given by the number of ways of 
forming such a loop, multiplied by a Boltzmann factor, 
e -8D/T therefore identify the entropy per loop as 

Sioop ~ — 2 log[4sinh(2 J/T)]. The apparent divergence 
of the perimeter (and therefore the entropy) at small J 
represents the breakdown of the ordered state which hap- 
pens at T c . The transition to the paramagnetic state 
occurs when the energy cost for adding two vertices to 
the loop (4D) is balanced by the entropy gain associated 
with adding an extra rectangular segment to the loop. 
This entropy gain is approximately Si oop ■ We may there- 
fore obtain an estimate of the transition temperature by 
setting AD/T C = S\ oop . The result is 



4sinh(2J/T c ) ~ exp(-2£>/T c ) 



(14) 



which differs from the exact result J2J only by the con- 
stant leading factor of 4. The real transition temperature 
is lower than that predicted by this method because in- 
teractions between the loops act to reduce the energies 
at large perimeters. 



Thus we interpret the transition to the paramagnet as 
deconfinement of localised composite excitations. The 
state becomes disordered when the loop size becomes 
comparable with the spacing between loops. 

The magnetisation and correlation lengths may also 
be calculated in similar series (see appendix El • The 
magnetisation, Mq , may be used to calculate the fraction 
of spins opposed to the mean spin, 



1-M 



-8D/T 



256 sinh 4 (2 J/T) 



(15) 



Assuming that the lowest-lying excitations are rectangu- 
lar domains with four excited plaquettes per domain, we 
expect a relation of the form n s = (n p / A) (n b / n p ) 2 . We 
see that this is true for small J/T [recall that (J/T) is a 
small number, although expansions about J = are not 
valid in general since we are in an ordered phase]. 

Thus we have arrived at the following picture of the 
thermodynamics of the SEV model. There is a density 
of excited plaquettes n p , which sit on the corners of over- 
lapping closed loops. The total loop perimeter is mea- 
sured in terms of the parameter nb: the point n b = (1/2) 
is the maximally disordered spin field, in which the ex- 
cited plaquettes are free. Away from the critical region 
(which is narrow), the excited plaquettes are nearly free 
in the paramagnetic phase: in the ferromagnetic phase 
then they are confined on corners of rectangular loops 
whose typical size is much smaller than their inverse den- 
sity. 

So far, our microscopic arguments have been purely 
thermodynamic: we have not considered any dynamics. 
In the next section, we investigate the dynamics of the 
SEV model. 



III. DYNAMICS 

This section contains the key results of this paper. We 
briefly describe the dynamics of the paramagnetic state, 
which are essentially independent of J. We then dis- 
cuss the onset of ordering as the temperature is low- 
ered through T c . We will find that supercooled states 
exist near T c , which are well described by a simple 'mo- 
bility field' picture for times shorter than their lifetime. 
We then discuss the extent to which these states can be 
regarded as metastable, and what determines their life- 
times. 

We begin with a very brief review of the dynamics 
in the paramagnetic state with T < T . Since J may 
be treated perturbatively in this regime, we write the 
Hamiltonian as in equation (|3J). Flipping a single spin 
involves flipping of the four plaquette variables adjacent 
to that spin. Thus spins adjacent to four unexcited pla- 
quettes have a flipping rate that is suppressed by a fac- 
tor e~ 8D / T . However, the flip rate of spins adjacent to 
exactly one excited plaquette are suppressed only by a 
factor e - 4D / T . The excited plaquettes mark mobile re- 
gions in which spin flips are rather likely. Thus the model 
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FIG. 3: Behaviour of the system after two quenches from in- 
finite temperature to just below and just above T c — 0.45D. 
The two plateaux correspond to the onset of activated dynam- 
ics, and equilibration in the disordered (glassy) state. How- 
ever, the disordered state is unstable below T c , as is clear at 
longer times. 



resembles kinetically constrained systems such as the FA 
modeljij. 

The relaxation time of the spins depends on the tem- 
perature and on the density of excited plaquettes, ac- 
cording to r ~ rip 1 e 4D / T . This arises from localised one- 
dimensional diffusion of pairs of excited plaquettes |2Ct 
l2l| . In equilibrium, we have n p ~ c = e~ 2 / J ', so the 
relaxation time diverges as c~ 3 . More precisely, we have 
the scaling relation 20J 

(aij(t w )aij(t + t t() )) c ,eq = f(c 3 t) (16) 

for the on-site autocorrelation function in equilibrium at 
a given value of c. This is strong glass behaviour in the 
classification of Angell 22] . 

We now turn to results for the SEV model below T c , 
where J may not be treated perturbatively. We dis- 
cuss the phenomenological similarities and differences be- 
tween this model and physical glass formers. We then 
interpret this behaviour with the aid of mean field rate 
equations. 



A. Existence of supercooled states 

We start this section by demonstrating how a super- 
cooled state may be formed below T c . We investigate the 
dynamics of the system by means of simulations which 
use a continuous time Monte Carlo algorithm 23J , with 
periodic boundary conditions. The number of excited 
plaquettes in each row and column is constrained to be 
even in this treatment, so the linear size of the system 
must be greater than (2/n p ) for reliable results. Some of 
the main features of the dynamics of the SEV model are 
shown in figure |3J We measure the energy density of the 



FIG. 4: Behaviour of the SEV model after a quench from 
Ti = 1.02T C to T 2 = 0.996T C with T c = 0.45D. We show 
concentrations of excited plaquettes, n„, and broken bonds, 
rib- The dotted lines mark ci,2 = e~ 2D ^ Tl - 2 . 



system with respect to the ground state: 

u = (H/N) +D + 2J (17) 

The dashed trace in figure shows the internal energy 
density after quenching to a temperature, Ti, such that 
T c < Ti < D. After an initial transient, the system 
ages in a power law fashion towards equilibration at 
u ~ 2De- 2D / Tl + 2J. The plateau at (u/D) ~ 0.2 is a 
characteristic feature of models with kinetic constraints 
(whether explicit 9] or emergent [lH Izjf). It represents 
the onset of 'activated' dynamics. The equilibration time 
for the paramagnet scales as a power of c = e~ 2D / T . All 
these features are seen in the model with J = 0, . 

In contrast, the full trace in figure [3] shows the be- 
haviour on quenching to a temperature T 2 satisfying 
T 2 < T c < D, but with T 2 close to T c . The behaviour 
resembles that of the quench to Ti, including apparent 
equilibration at u ~ 2De~ 2D / T2 + 2 J. However, this is 
state is not stable, and the energy falls further at longer 
times. This behaviour resembles that of glass-formers, 
where the state on the lower plateau would be called a 
supercooled liquid. The behaviour is also qualitatively 
similar to that observed in [l2|. 

In order to focus on the supercooled states, we show 
further simulations in figure 0] The system is cooled 
through T c in a single small step. We plot n p and rib 
as a function of time. From the plot of n p , we see that 
the density of free excitations responds relatively quickly 
to the change in temperature: it falls from c\ = e -' 2D / T i 
to C2 = e~ 2D / T2 , where it appears to stabilise. 

The response of rib to the cooling is much slower. Re- 
call that this correlation function measures the clustering 
of excited plaquettes. This clustering is a much slower 
process than the creation and annihilation steps leading 
to a change in the concentration of excitations. Look- 
ing at the late times in figure 0] when the clustering does 
start to occur, the system ages towards the ferromagnetic 
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FIG. 5: (Top) fraction of broken bonds against waiting time in 
the supercooled state: T = 0.996T C , T c = 0.444D. (Bottom) 
Autocorrelation function C(t) — {<Jij{t + t w )uij(t w )} for the 
waiting times marked in the top figure. The two-time func- 
tion appears stationary until the system starts to nucleate at 
around 3 x 10 7 Monte Carlo sweeps. 



state with both n p and rib falling together. Taking the 
two traces in figure U together, we see that there are two 
separate timescales: one is associated with changes in n p ; 
the other with changes in rib- 

Turning to the supercooled state itself, it is clear from 
figure that it has n p ~ c = er 2D l T and rib ~ (1/2). 
This resembles closely the state that would be formed if 
we set J = in the Hamiltonian. Thus the effect of the 
interactions between plaquettes (the term proportional to 
J in the Hamiltonian) is to set the lifetime of the super- 
cooled state. The properties of the state itself are inde- 
pendent of J. We conclude that the supercooled state in 
the SEV model can be well described by the much simpler 
plaquette model of equation J3J : a kinetically frustrated 
model with trivial thermodynamic properties. This is the 
assumption made when describing glass-formers by sim- 
ple models of dynamical heterogeneity 0, 0|: in the SEV 
model, this assumption seems to be reasonable. 

A key property of supercooled states is that two-time 
correlation functions should be stationary within the su- 
percooled state. That is, expectation values of the form 
(A(t w )A(t + t w )) should be independent of t w as long as 
(t + t w ) is less than the lifetime of the states. In fig- 
ure El we show that the supercooled state with n p = c 
and rib = (1/2) na s this property. Since the excited pla- 
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FIG. 6: Nucleation times at T c = 0.50D and T c = 0.44D, 
measured by averaging over several runs similar to those in 
figure [5] The temperature increases from left to right, but 
note that the zero of the logarithmic scale is at T c , away to 
the right. The dashed lines are fits to the form of equation 
11811 using the same value of 7 = 0.17. This form captures the 
qualitative features of the curve: a power law increase near 
T c , and an exponential divergence at low temperatures. 



quettes are uncorrelated in this state, it may be prepared 
manually, without the need for simulation. Further, the 
single spin autocorrelation function in the supercooled 
state is the same as that in the model with J = at 
the same temperature. It therefore obeys the dynamical 
scaling law of equation (|16fl . and is independent of J. 

As a final comment on figure [3 we note that the cri- 
terion that the supercooled state should be long-lived 
is fulfilled, since stationarity holds over timescales much 
longer than the single-spin relaxation time. Thus the 
data in that figure is consistent with the picture of 
a supercooled state that appears to equilibrate in a 
metastable basin. 



B. Lifetime of supercooled states 

Having identified a supercooled state at temperatures 
near T c , we now proceed to discuss its lifetime. We pre- 
pared states with n p = c and rib = (1/2) and measured 
the time taken for rib to fall to 0.45. Results are shown in 
figure El The lifetime gets very large both near T Cl and 
in the limit T — ► 0. 

Near T c , the states are supercooled. At very low tem- 
peratures, the lifetimes are similarly long, but in this 
case they are of the same order as the spin relaxation 
time. These states are not equilibrated in a metastable 
basin; rather, their long lifetimes reflect the drastic slow- 
ing down of all timescales as the temperature is reduced. 

From figure El we conjecture that the nucleation time 
has the form 



7e 



2D/T a 



1 - (T/T c ) 



(18) 
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where 7 is a (microscopic) rate that depends only weakly 
on T and T c . Physically, the nucleation rate is suppressed 
by a Boltzmann factor that arises from the activated dy- 
namics, and by a factor of A/x/T ~ (T c /T) — 1, where 
A/i is the free energy difference between ordered and dis- 
ordered phases. 

In other words, behaviour near T c is characterised by 
a separation of the nucleation time from the relaxation 
time, r. The minimum in the nucleation time occurs at 

T x = T c [l - (TJAD) + 0{T c /Df] (19) 

The relaxation time in the paramagnet given ap- 
proximately by r = 0.14e 6I) / T , so at the mini- 
mum we have r* ~ 0.14e (6I)/T <= )+(3/2) and r* uc ~ 
0.17(4:D/T c )e (6D / T ^ +1 . The result is that r* uc > r* as 
long as T c < T Q ~ D, which is the regime of interest in 
this section. The two timescales are well separated for 
all temperatures between T x and T c . 

This separation results from the small amount of free 
energy that is released on ordering. States in which these 
times are well-separated are 'supercooled', in the sense 
that the degrees of freedom associated with the relaxation 
time appear to equilibrate in a state that is known to be 
unstable at long times. 

At very low temperatures, we see that r nuc will become 
smaller than the extrapolated relaxation time tj = o ~ 
e 6D/T _ rp ne resu j^ ihax the physical relaxation time at 
low temperatures is smaller than e 6D / T . The dynamics 
of the aging state are faster than those of a similar state 
with J = 0. 

These results may be interpreted in the picture of the 
model as a combination of a zero temperature dynami- 
cal fixed point and a finite temperature thermodynamic 
singularity. The spin relaxation time is controlled by 
the activated dynamics associated with the zero temper- 
ature fixed point. It is large compared to microscopic 
timescales, but small compared to the lifetime of the 
supercooled state. That lifetime is very long near the 
thermodynamic transition: the slowing down is due to 
the small free energy difference between paramagnetic 
and ferromagnetic states. We comment here that 'soft 
modes' at large lengthscales are not relevant to the be- 
haviour observed in simulations, due to the narrowness 
of the critical region. 

While supercooled states are familiar in systems with 
first order phase transitions, they are not usually ob- 
served near second order transitions, such as the one 
discussed in this work. In first order systems, the nu- 
cleation time may be predicted by thermodynamic argu- 
ments. The free energy of a (d-dimensional) droplet of or- 
dered phase in a disordered background is approximately 
ar 1 — A/j,r d where r is the linear size of the droplet, a 
the surface tension and A^ is the free energy difference 
between the two phases. Thus, nucleation requires the 
formation of a droplet of linear size r* ~ er/A/i, with 
an associated free energy barrier proportional to A/x -1 
(in two dimensions). The nucleation rate is therefore 



suppressed by a factor e~ a / TA ^. This is exponential 
suppression of nucleation. 

From (|18fl . we see that the SEV model has linear 
suppression: the nucleation rate is proportional to A/i. 
This second order system has weaker suppression than 
that predicted for first order systems. Since the transi- 
tion in the SEV model is second order, there are pro- 
cesses by which the bulk of the system may be continu- 
ously changed from paramagnet to ferromagnet, without 
a large free energy barrier. These processes are slow be- 
cause they require co-operative motion of many spins, 
but the exponential slowdown that would result from a 
high energy intermediate state is not observed. While the 
phenomenology of the SEV model resembles that of first 
order systems, the lifetimes of supercooled states tend to 
be shorter than those in systems with diverging free en- 
ergy barriers near T c . In this respect, the SEV model is 
an imperfect model of a glass-former. However, we argue 
that the free energy barrier between ordered and disor- 
dered states in first order systems should mean that the 
supercooled states are less affected by the critical point 
than those of the SEV model. Thus, if the thermody- 
namic singularity is largely irrelevant in this model, then 
we expect it to be even less relevant in similar first order 
systems. 



C. Rate equation approach and aging behaviour 

In order to understand the results of the previous sec- 
tion, we give a brief discussion of the aging behaviour of 
the system. We parameterise this behaviour in terms of 
mean- field rate equations for the observables n p and n^. 
This will provide further evidence that the supercooled 
states are characterised by fast dynamics for the concen- 
tration of excited plaquettes, n p , combined with much 
slower dynamics for their spatial ordering (measured by 
rib). 

In paramagnetic states with T < D, the equilibrium 
value of rib is 0.5 and that of n p is c ~ e~ 2D l T . Aging 
towards equilibrium occurs at rib = 0.5, with n p ~ t~ 0A5 
(a sample trace is shown in figure but the observed 
exponent is independent of temperature, as long we work 
between T c and T Q ). The rate is limited by the slow 
diffusion of excited plaquettes (there is no simple diffusive 
process for isolated excitations). 

There are two regimes for the aging towards the fer- 
romagnetic state. We define T x as the temperature at 
which the nucleation time is minimal (recall figure 
As shown in figure [21 after quenching to T near T c 
(T > T x ), the system appears to equilibrate at n p = c 
(with rib = 0.5), before aging towards the ferromagnetic 
state. In figure [S] we show a similar scenario, but con- 
tinuing to slightly longer times. At these long times, the 
system ages with the ratio (n p /2rib) ~ c. The decrease 
of rib happens on the long timescale given by r nuc . How- 
ever, the plaquette concentration, n p , is reacting on a 
much faster timescale. The interpretation is that n p has 
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FIG. 7: Plot showing aging of the paramagnet. We plot n p 
as a function of time after a quench from infinite temperature 
to T — 0.286D, with J = 0. The dashed line is a guide to the 
eye, corresponding to a decay proportional to i -0 ' 45 . 
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FIG. 9: Aging at T = 0.364D, with T c = 0.444D, so that 
T < T x . The ratio n p /2nt deviates from n p at a value around 
e -2D/T a ^ go ^ ne S y S t em does not reach the supercooled state 
with n p = c = e' 2D/T . 
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FIG. 8: Aging of the system at a temperature T = 0.995T C . 
We have T c = 0.451D so that T > T x . We plot n p and n p /2n b 
after a quench from infinite temperature. Both quantities 
are equal until the supercooled state becomes unstable. At 
these long times, n p decreases, but the ratio n p /2nt remains 
constant with a value approximately equal to c = e~ 2D ^ T . 



stabilised in an environment set by the current value of 
n b . It seems that the stable value of n p is around 2cn b . 

If instead we quench to a temperature below T x , there 
is no apparent equilibration at n p ~ c. Instead, the ef- 
fect of finite J becomes apparent when n p ~ e -2£>/T c > 
c. There is a transient effect as rib begins to change, 
but the long time behaviour again has a constant ratio 
(n p /2n b ) — 0.8e~ 2D / Tc (see figure Simulations indi- 
cate that the constant value of {n p /2n b ) seems to depend 
on T c as shown, and to vary only weakly with T. Further, 
the exponent associated with the aging appears to be sim- 
ilar to that in the paramagnetic state (around 0.45). The 
natural interpretation is that the aging of these states is 
controlled in the same way as the aging of the param- 
agnet: by the slow diffusion of excited plaquettes. The 
dynamics of rib are now faster than those of n p , and the 



spatial ordering of the excited plaquettes is responding 
faster than the plaquette density. 

This description of the aging behaviour is consistent 
with the following conjectured rate equations 

e iD ' T d t n p = X p n 2 p 2 (2cn b -n p ) (20) 
e 4D/T e 2D/T cdtnb = Xb n b [l-(T/T c )]{X b -2n b ) 

(21) 

where X b = min(l, n p Ae 2D ^ Tc ). The significance of this 
quantity is that n b would equilibrate to X b /2 if the con- 
centration of plaquettes were constrained to be equal to 
rip. The exponent 2.2 in equation l|20|) is fixed by the 
exponent 0.45 for the power law decay of the energy in 
the paramagnetic phase (as shown in figure 01 ■ The the 
overall scaling of time with temperature in that equation 
is also fixed by the aging of the paramagnetic states. The 
scaling of time in l|21|) is determined by the scaling of the 
nucleation time in equation (|18f) . 

The adjustable parameters of the theory are there- 
fore the rates A p and Ah, and the constant, A. The 
two rates are microscopic frequencies reflecting the co- 
operative motion of the spins required to change n p or 
n b . The constant A determines the ratio of rib and n p 
when the aging is controlled by plaquette diffusion (as in 
figure EJ). The instability of the supercooled state means 
that Ae 2D / Tc < e 2D / T : the data of figure is consistent 
with (1/A) ~ 0.8, as mentioned above. 

We make no attempt to justify these rate equations 
on microscopic grounds. For example, the exponent 2.2 
in equation l|20() is a non-trivial decay exponent for an 
annihilation-diffusion problem in which diffusion of sin- 
gle excited plaquettes is suppressed by the kinetic con- 
straint. We imagine 'integrating out' all the microscopic 
degrees of freedom: the effect of the complicated fluc- 
tuation effects appears only in this exponent. However, 
while their microscopic origin is unclear, the interesting 
feature of these equations are (1) the different tempera- 
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ture scaling of the times associated with the two degrees 
of freedom, and (2), the presence of points at which one 
degree of freedom is not changing. The first feature leads 
to a separation of timescales in the problem. When this is 
combined with the second feature, the apparent metasta- 
bility of the supercooled states becomes possible. In this 
case, the faster degree of freedom is n p , which equili- 
brates at 2crib on a timescale that is fast enough that rib 
can be considered to be constant. The aging of the su- 
percooled state then has a timescale determined by the 
rate equation for rib- That degree of freedom is trying 
to reach apparent equilibration at n p Ae 2D ^ Tc /2, but n p 
is the faster degree of freedom, and is moving the target 
downwards just as fast as n& decreases towards it. This 
results in the aging at constant rib/rip that is shown in 
figure |SJ An exactly analogous process is taking place in 
figure^! except that n p is now the slow degree of freedom. 

The central point of the above argument is that the 
slow degrees of freedom set a 'target value' for the fast 
ones, at which the fast degrees of freedom appear to equi- 
librate when viewed on the fast timescale. This is the 
sense in which the states discussed in the previous sec- 
tion are 'supercooled'. Their lifetime is then set by the 
slow degrees of freedom, and this lifetime is much longer 
than relaxation timescales in these states. To understand 
the correlations in the aging state, it would be desirable 
to study the thermodynamics of the excited plaquettes 
while working at fixed rib- However, this is well beyond 
the scope of this paper. 

Note that there is no provision in equations (|2(J[I and 
(I21|l for equilibration in a ferromagnetic state. However, 
the large internal energy difference between paramag- 
netic and ferromagnetic states means that this equilibra- 
tion is not observed on timescales that are accessible in 
simulation. Therefore, equations l|20(l and l|21|) appear to 
be valid over timescales that are several orders of magni- 
tude longer than the lifetimes of the supercooled states. 

It is simple to verify that the qualitative behaviour of 
figuresEHand l7l9l mav be fitted by equations lj2"U|) and (|2"T|> . 
with appropriate values of (A p , At, A). See figure HUI in 
which we show reasonable agreement. Note however, that 
the onset of nucleation from the supercooled state is more 
sudden than that predicted by the mean field equations. 
The initial ordering is slower than the power law pre- 
dicted by these simple rate equations. 

The requirement of different parameters to fit the dif- 
ferent simulations in figure ITUl show the rather simplistic 
assumptions for the temperature scaling in the mean field 
equations. That is, the linear suppression of the nucle- 
ation rate with (T c —T) is valid only near T c , necessitating 
adjustment of \b at smaller temperatures. We have al- 
ready commented that A will be temperature dependent, 
but that this dependence has a weaker effect on Xb than 
the exponential dependence of that quantity on T c . 

The fits of figure ^3 using these mean field equations 
support our interpretation of figure El in terms of purely 
dynamical effects that do not depend on thermodynamic 
quantities like free energy barriers and spinodal points. 
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FIG. 10: Fit to the data of figures H (+ symbols) and HO (x 
symbols) using mean field rate equations 121)1211 . Note that 
the left and right panels each show data from two runs, show- 
ing how the different regimes are captured by the same rate 
equations. Parameters are (A p , \t, A) = (350, 17, 0.85) for the 
data of figure H and (X p , X b ,A) = (350, 17, 1.20) for that of 
figure |U] 



IV. DISCUSSION AND CONCLUSIONS 

We have shown that the SEV model can be interpreted 
in terms of a high temperature state in which excitations 
are free and point-like, and a low temperature phase in 
which these free excitations are confined into composite 
objects, with a characteristic size that is smaller than 
their spacing. The dynamics of the paramagnet are well 
described in a 'mobility field' picture similar to the FA 
model. 

Near the transition to the ordered phase, 'supercooled' 
paramagnetic states have long lifetimes. Two-time cor- 
relation functions show stationarity in these states, and 
both their dynamics and their thermodynamics are well- 
described by the plaquette model. However, these states 
have finite lifetimes, beyond which stationarity is lost 
and it becomes clear that they are unstable to the fer- 
romagnetic state. This situation resembles the situation 
in physical glass-formers, except that the presence of a 
first order transition in those systems means that the life- 
times of the supercooled states diverge much faster near 
T c than in the SEV model. 

We have also shown that the mean field rate equations, 
(|20|l and (|21|l , are a suitable framework for describing the 
out of equilibrium (aging) behaviour of the SEV model. 

We end this work with some comments about the sig- 
nificance of these results in the context of the literature. 
Comparing with the work of Cavagna et al. fl^ , we note 
the similarity between the phenomenology of their model 
and the SEV model. However, the exact solution of the 
eight vertex model, and the understanding of the dynam- 
ics of the paramagnetic state that arises from previous 
work on the plaquette model allow us to investigate the 
behaviour from a more microscopic viewpoint. 

We illustrate this with three points. Firstly, the in- 
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ternal energy of the SEV model changes very rapidly 
around T c . From simulation evidence alone, we might 
(erroneously) conclude that the transition was first or- 
der. However, we know from theoretical considerations 
that the transition is second order. This knowledge is 
important when discussing the possible metastability of 
supercooled states. Secondly, the power law suppression 
of nucleation near T c in the SEV model results from the 
fact that our transition is second order. For a first order 
transition we would expect an exponential suppression. 
The nature of this suppression in [12 ] does not seem to be 
clear. And finally, we are able to identify the minimum 
in figure as arising not from the crossing of a spinodal, 
but rather from a crossover in timescales. 

We would also like to point out some similarities be- 
tween the paramagnetic phase of the SEV model and the 
models which are conjectured to be controlled by the be- 
haviour of entropic droplets 0)113 • As mentioned above, 
the invariance of the Hamiltonian at J = under flipping 
a whole row or column of spins leads to a large (but non- 
extensive) ground state entropy. We have shown that the 
introduction of J is largely irrelevant above T c , and in the 
supercooled states. Therefore, we may interpret the low 
temperature paramagnetic states (T c < T <C T a ) as a 
'mosaic' of the many J = ground states. The 'droplets' 
are referred to as 'entropic' in this scenario. This name 
arises because there does not seem to be an energetic ar- 
gument for their stability, so it assumed that they are 
stabilised by some entropic mechanism. In the plaque- 
tte model, the borders between 'droplets' of each ground 
state are not one dimensional as one might expect, but 
rather arise from point-like excited plaquettes. This situ- 
ation was alluded to in the recent paper of Bouchaud and 
Biroli 25]. Since there is no perimeter for a surface ten- 
sion to couple to, the entropic forces are strong enough 
to stabilise the 'mosaic' state. This is the situation in the 
SEV model at finite J, but T > T c . 

If we accept the plaquette model as a realisation of 
the entropic droplet picture, it is interesting to note that 
the only fixed point in that model is at zero temperature. 
That is, although the entropy falls rapidly as the temper- 
ature is reduced through T Q , any extrapolation that leads 
to S = at a finite temperature is not valid. Rather, the 
entropy remains regular at all temperatures (even those 
in which the glassy phase is completely unstable to the 
crystal). In fact, iS g i aS s > ^crystal a H the way down to 
T = (where S'giass is the entropy associated with the 
plaquette model and S C rystal is the entropy of the SEV 
model with the same value of D and at the same tem- 
perature). In this scenario the Kauzmann paradox |22| 
is seen as arising from an unphysical linear extrapolation 
of the glassy entropy. 

To conclude, we have shown that the plaquette model 
limit of the SEV model describes its behaviour in both 
stable and supercooled paramagnetic states. This model 
resembles both the mobility field description of glassy 
systems (as exemplified by the FA model) and the en- 
tropic droplet picture. However, there is no finite temper- 



ature fixed point in the theory of the glassy states; thus 
Kauzmann's paradox is avoided. Taken together, these 
results are further evidence that theories without thermo- 
dynamic singularities at finite temperature are suitable 
for describing glassy states 0, 0, E| . 



Acknowledgments 

We thank A. Cavagna, D. Chandler and S. Whitelam 
for discussions. We acknowledge financial support from 
EPSRC Grants No. GR/R83712/01 and GR/S54074/01, 
University of Nottingham Grant No. FEF 3024. 



APPENDIX A: LOW AND HIGH 
TEMPERATURE EXPANSIONS OF THE EIGHT 
VERTEX MODEL 

In this appendix, we consider the series for the free 
energy given by Baxter [11^: 



(//T) = -logc / 



yv x m ((q/x 2 ) 



2\ro 



l) 2 (l-X m ) 2 



m(l - q 2m )(l + x m / 2 ) 



( A1 ) 

where q, x 2 and d depend on the model parameters 
(D/T) and (J/T). The parameter z in Baxter's calcu- 
lation is 1 in the SEV model since the Hamiltonian is 
invariant under 90° rotations. The dependence of q and 
x 2 on the original parameters (D, J) is rather indirect: 
the main task of this appendix will be to derive simple 
relations between (q, x 2 ) and (D/T, J/T) in the ferro- 
magnetic phase. 

However, we first consider the contribution of the logc' 
term to the free energy. In the paramagnetic phase we 
have 



(/pm/T) ~ logc' = log [e D / T cosh 



(A2) 



which is the result quoted as 10, and may be used to 
calculate properties of the paramagnetic phase. How- 
ever, in the ferromagnetic phase we have d — e^ D+2J ^ T , 
so taking only the leading term leads to the trivial result: 
Tip = Tib = 0. There are no excitations in the ferromag- 
netic phase at this order. We must therefore estimate the 
parameters x and q in this phase. 

The prescription for calculating these parameters is 
given by Baxter [19j, but we give a brief review. The 
ratios D/T and J/T are used to calculate four param- 
eters (a, b, c, d). The partition function is symmetric 
under interchange of the four quantities (a ± b, c± d). 
We may therefore map the parameters into the principal 
regime, to give (a', b' , d, d') satisfying d + d' > d — d' > 
a' + b' > a' — b' > 0. These parameters are used to cal- 
culate A = (-(d 2 + d 12 ) + (a' 2 + b l2 ))/{2dd' + 2a'b') and 
1 =(dd'/a'b'). 

In the ferromagnetic phase we have 7 = exp(4D/T). 
We work at T c < D so 7 is a very large number, and we 
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arrive at 



A = cosh(4 J/T) [1-^(1 



1 



: i+e>( 7 - 2 )] 

(A3) 



cosh(4 J/T) ' 

The next step is to calculate the parameter k, from: 

7 (l + fc 2 ) = A 2 (l+ 7 ) 2 -(l+ 7 2 ) (A4) 
and the result is 

-4D/T 

k = = + 0( 7 " 2 ) (A5) 

2 sinh 2 (2 J/T) [l + cosh(2 J/T)] w ' V ' 

We have < k < 1, and fc is the elliptic modulus with 
nome g. That is, 



logg 



dt 



dt 



/o vT^Vi _ p^i 7 y/TT¥ Vk 2 + t 2 

(A6) 

For the series of IjAlfl to converge quickly, we require 
q < 1: in that case (|A6I) reduces to q ~ (fc 2 /16), and 
therefore we have the approximate relation 



3 -8D/T 



128 sinh 4 (2 J/T) 



(A7) 



whose condition for validity is that k be small, which 
requires that sinh(2 J/T) < e~ 2D/T . This is the condi- 
tion that we are well inside the ferromagnetic phase, as 
is clear from @. 



In order to evaluate the terms in (|A1|) , we also require 
an approximate form for x. The definition of x is 



dt 



irdt 



o VT^ tVI - Pi 2 V^VT+tVPTt 5 

from which the relation a; 2 < g is clear. This integral 
may be expanded in a series around (l/\/jk) — 0. Equa- 
tion 1A5|I shows that this parameter is small if J <C T. 
The result is that 



(q/x 2 ) - 1 ~ 4 sinh(2 J/T) 



(A9) 



Substituting into (|Aip> . and ignoring all terms with m > 
2, we arrive at 



P -8D/T 

/ fm ~ —(D/T) (J/T) -r^r-TT^ 



16 sinh (2 J/T) 



(A10) 



which gives the result l|ll|) , qualified by the validity con- 
dition 



' 2D/T < sinh(2 J/T) < 1 



(All) 



from which we note that this is not an expansion about 
T = 0, but rather a useful approximation to the free 
energy in this part of the parameter space. 
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